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Abstract. We develop and implement an algorithm for the quantitative 
characterization of cluster dynamics occurring on cellular automata defined on 
an arbitrary structure. As a prototype for such systems we focus on the Ising 
model on a finite Sierpsinski Gasket, which is known to possess a complex 
thermodynamic behavior. Our algorithm requires the projection of evolving 
configurations into an appropriate partition space, where an information-based 
metrics (Rohlin distance) can be naturally defined and worked out in order to 
detect the changing and the stable components of clusters. The analysis highlights 
the existence of different temperature regimes according to the size and the rate of 
change of clusters. Such regimes are, in turn, related to the correlation length and 
the emerging "critical" fluctuations, in agreement with previous thermodynamic 
analysis, hence providing a non-trivial geometric description of the peculiar 
critical-like behavior exhibited by the system. Moreover, at high temperatures, we 
highlight the existence of different time scales controlling the evolution towards 
chaos. 
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1. Introduction 

The effect of inhomogeneity on the criticaf behaviour of magnetic systems has been 
considered in various contexts (e.g. disorder, coupling randomness, quasiperiodic 
structures); in particular, discrete-spin models defined on fractal topologies possess 
critical properties significantly different and richer than those found for translationally 
invariant systems [3 HJ [3J U E3 ■ 

The interest in fractal structures is not purely theoretical: many condensed- 
matter systems display strong nonuniformity on all length scales and can therefore be 
characterized as fractal objects; examples include the backbone of percolation clusters, 
aggregates obtained from diffusion-limited growth processes, and absorbent surfaces. 

One of the most known fractals is the Sierpinski gasket (SG), which, due 
to its exact decimability, allows analytical approaches; in particular, by means of 
renormalization group techniques, it was proved that the Ising model on the SG 
exhibits phase transition only at zero temperature, while at any finite temperature the 
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system breaks into domains and loses long-range order [TJ [3J [3J [7] . While this result 
was found in the thermodynamic limit, at the mesoscopic sizes peculiar and interesting 
thermodynamic properties arise [8- • More precisely, the Ising model defined on a finite 
SG exhibits critical-like features at nonzero, low temperatures and the solution found 
in the thermodynamic limit turns out to be a poor approximation. This anomalous 
behavior has been investigated from a thermodynamic point of view and derives from 
long-range, slowly decaying correlations at low temperatures [9]. 

The way such a behavior is reflected by the evolution of spin configurations 
is an item so far overlooked (even if, for rectangular lattices, the idea of studying 
cluster dynamics may be traced back to Peierls and Griffith [TQl [11] ). Given the 
importance of the dynamics of Ising-like clusters in many research areas, from 
condensed matter to biological systems Q21 EH H] , the definition and the development 
of proper tools for this kind of analysis would be very useful. Moreover, it would be 
particularly intriguing for inhomogeneous substrates, due to the emerging non-trivial 
thermodynamic behavior; on the other hand, it is just on such structures that the 
definition of a proper metrization or evolutionary dynamics can be more awkward. 

The interest in cluster mobility actually extends to an extremely wide class of 
models: indeed the Ising model on the SG may be seen as a particular realization of 
cellular automata on graphs, i.e. discrete time dynamical systems assigning to each 
node of a graph G a value chosen in a alphabet K, along a rule depending only on a 
finite neighborhood at previous time [15) . 

In this work we aim to introduce and develop proper algorithms for the 
quantitative characterization of cluster dynamics on graphs, and we use this approach 
for the Ising model on a finite SG, meant as a prototype of cellular automata 
on graphs. The procedure requires a projection of evolving configurations into an 
appropriate partition space, where an information-based metrics (Rohlin distance) 
and a method measuring the effective emergence of configurational novelty (reduction 
process and amplification parameter) may be naturally defined and worked out in order 
to focus the changing and the stable components of configurations. The algorithmic 
implementations of Rohlin distance and related quantities are deeply affected by the 
topological features of the substrate (T5J- For instance, in previous implementations 
designed specifically for automata on regular lattices, the very passage from one to two 
dimensions yields a much higher order of complexity, and could not be exported on 
different substrates [T71 [HI HH [20] . On the contrary, the algorithm developed here can 
be directly applied to generic cellular automata for which the metric characterization 
of cluster dynamics gets feasible; a brief description is given in the appendix B. 

Our investigations on the Ising SG highlight the existence of two "critical" 
temperatures, Tj and Tu, demarcating three main regimes which recover, both 
qualitatively and quantitatively, the results of the previous thermodynamic analysis; 
in addition, within the above mentioned regimes, we obtain a more detailed dynamical 
and geometrical characterization. In particular, at very low temperatures (T < Tj), a 
long-range order is established, and the few small-sized clusters display poor overlap 
from one time step to the next one: this makes the distance close to zero, and the 
reduction ineffective. As the spot sizes start to increase, the slowness of the evolution 
is still such that both non-similarity and overlap between successive configurations 
rapidly grow. At greater temperatures (Tj < T < Tu), large-scale correlations start 
to decay, clusters of all sizes appear, and overlap gets easier; then, for T > Tu, any 
trace of order has vanished and, even if overlaps are very important, the complexity 
of magnetic pattern is irreducible. The progression of such different kinds of non- 
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similarity and dynamical overlap is well described by our parameters (distance, 
amplification and intersection, see below). We will also evidence the existence of 
different scales controlling the disappearance of local and correlated order. Thus, the 
phenomenology provided by our method, confirms with a deeper geometric insight the 
peculiar critical-like behavior exhibited by the system 

We finally notice that, in view of future extensions to non-equilibrium situations, 
we will adopt a microcanonical dynamics |21) . which allows implementations even in 
the presence of temperature gradients. 

In the next two sections, after recalling basic notations on graphs, we review some 
facts concerning the thermodynamic properties of the Ising model on the Sierpinski 
gasket (Sec. [2]), and the microcanonical dynamics working as evolutionary dynamics 
(Sec. [J]). Then, we introduce general procedures for cluster identification and reduction 
(Sec. and we show our results on the SG (Sec. [5]). Finally, we present our conclusions 
and perspectives (Sec. [5]). Technical remarks about partitions on graphs can be found 
in the appendices. 

2. Ising model on the Sierpinski gasket 

A generic graph G is mathematically specified by the pair {A, T} consisting of a non- 
empty, countable set of points A joined pairwise by a set of links T. The cardinality 
|A| = N of A represents the number of sites making up the graph, i.e. its volume. 
From an algebraic point of view, a graph is completely described by its adjacency 
matrix A. Every entry of this off-diagonal, symmetric matrix corresponds to a pair 
of sites, and it equals 1 if and only if this couple is joined by a link, otherwise it is 0. 

Here we consider the Sierpinski gasket which can be built recursively with the 
following procedure: the initial state (Go) is a triangle and the gth stage G g is obtained 
joining two of the three external corners of three G 5 _i to form a bigger triangle (see 
Fig. QJ. In this way, the volume of G g is 3(3 9 — l)/2. The Gasket is obtained as the 
limit for g — > oo of this procedure. 

2.1. Thermodynamic properties of the Ising model on the SG 

The Ising model on a generic graph G = {A, T} is defined associating the spin variable 
Si — ±1 to every site i of the graph, and considering a nearest-neighbours interaction 
between points i and j, such that Aij = 1. The Hamiltonian is therefore 



where s = {s^jigA denotes the magnetic configuration of the system and the coupling 
J is assumed to be the same for any couple; in the following we will set J = 1. 

As it is well known, a magnetic model defined on a finite lattice cannot exhibit 
critical behaviour at nonzero temperatures; critical features can only emerge when the 
underlying lattice becomes infinitely large, i.e. in the thermodynamic limit. Another 
necessary condition in order to have a nonzero critical temperature concerns the 
topology of the (infinite) substrate: for Euclidean structures it was rigorously shown 
that the dimension must be larger than 1. Analogously, it has been shown that the 
discrete symmetry Ising model on finitely ramified fractals cannot have a nonzero 
critical temperature p3 [3J However, the critical behaviour of one-dimensional 

systems and finitely ramified fractals can be markedly different, since for the latter it 
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Figure 1. (Color on line) Sierpinski gasket of generation 2 (a), 3 (b) and 4 (c), 
with volume N = 6, N = 15 and N = 42, respectively. 



is further governed by additional geometric aspects such as ramification, lacunarity, 
and connectivity pQ. 

The thermodynamic properties of the Ising model on the SG were studied in 
details in p] [21 [3l [H HI El [23] , where it was shown that its scale-invariant, fractal 
structure leads to highly cooperative correlations and, at sufficiently low temperatures, 
the correlation length £ becomes extremely large and slowly decaying (compared with 
the one-dimensional case), so that any system with size smaller than £ displays long- 
range order. Indeed, one can define an apparent magnetic transition temperature T\ as 
the point where (to 2 ) = 0.5; for a system of generation g (g ^> 1), T\ ~ AJ/ ln(4g), for 
example, for generation g — 5 and g — 6 one finds T\ « 1.335 and T\ ~ 1.259, 
respectively. Hence, as the system size is enlarged, T\ diminishes slowly; more 
precisely, being TV = 3(3 5 + l)/2 the number of nodes for a gasket of generation 
g, one has T x ~ 1/ ln[ln(7V)]. 

On the other hand, the specific heat c does not display any anomaly associated 
with long-range order, yet it exhibits a peak at a temperature T 2 « 2 J [8]. Conversely, 
the "reduced" specific heat ~ cT 2 , (basically the derivative of the magnetic energy 
per link with respect to the inverse temperature [5]), evidences a qualitative difference 
between the SG and a one-dimensional system, since for the former it exhibits a peak, 
while for latter it grows monotonically as the temperature is increased. Indeed, for 
the SG, as T is increased from small values, the energy decreases rather slowly as a 
result of a relatively large cost due to the large (bulk) coordination number; at larger 
temperatures it becomes progressively easier to reduce the energy and this reflects 
the fact that the large fluctuations begin to develop; finally, since the energy must 
ultimately vanish, the rate of change gets smaller and smaller [9] . 

In conclusion, the SG displays a non-trivial thermodynamic behavior which can be 
summarized as follows. At T < T\ a long-range order is established; for T\ < T < T 2 a 
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short-range order is still present with large fluctuations on all length scales. At larger 
temperatures T > T 2 any trace of order has disappeared and a paramagnetic state is 
approached. 

3. Microcanonical Dynamics 

When studying transport properties, microcanonical dynamics are usually chosen as 
they allow to describe an isolated system, or its isolated bulk, without any assumption 
on the equilibrium state between the system and the surrounding. 

Here we adopt a recently introduced microcanonical dynamics [21] , which features 
a high degree of flexibility, being ergodic in any temperature range and implcmcntable 
on a generic structure, even in the presence of disorder. On regular lattices (e.g. 
cylinder, torus) such a dynamics has already been shown to be able to lead the system 
to thermalized states compatible with those expected from a canonical dynamics and 
to allow the study of out-of-equilibrium properties [3TJ [51] . 

Although we will focus only on equilibrium regimes, the reason for choosing 
this dynamics is twofold: First we test its reliability on an inhomogeneous structure; 
second, we pave the ground for the study of transport properties on such a substrate. 

In the following, we briefly resume how it works having in mind as substrate a 
generic graph G described by the adjacency matrix A. 

For each pair of connected sites, namely each link such that Aij — 1, besides 
the magnetic energy E™ = JsiSj, we introduce a local kinetic energy Eij > which 
is, in principle, unbounded. Now, the dynamical rule proceeds as follows: 

(i) Start from a (discrete) distribution of energies {Sjj}ij=i jv; 

(ii) choose randomly a link i ~ j; 

(hi) extract one over the possible four spin-configurations for the couple of sites i, j, 
and evaluate the magnetic energy variation AE 7 ^ induced by the move; 

(iv) if AEfj < 0, accept the move and increase the link energy £Jy of AE^. When 
AE"j > 0, the move is accepted only if E^ > AE™ and the link energy is 
consequently decreased of AE™ ; otherwise the move is not accepted and the link 
energy is not updated. 

It is worth remarking that AE™ allows for energy variations occurred on the link 
i ~ j as well as on those pertaining to links adjacent to sites i or j: 



where Si — 1 if the ith spin has undergone a spin flip, otherwise it is zero. It is therefore 
clear that, due to the discreteness of the system and to the fact that J is constant 
over all links, both E*™ and E^ are discrete variables: the former can only assume 
two different values corresponding to the aligned and non aligned configurations of the 
adjacent spins i and j; the latter can only assume integer values deriving by proper 
combinations of the pertaining AE™. 

We also notice that Eij works as an additional degree of freedom and the above 
dynamics conserves the total energy given by the following Hamiltonian function 
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As shown in |21) . the magnetic and kinetic energies result to be non-correlated: 
this allows a natural definition of temperature at equilibrium, which, in a very natural 
way, depends only on the average kinetic energy. In fact, the link energy satisfies the 
Boltzmann distribution exp(—(3Eij), and the fitted constant j3 — 1/T just corresponds 
to the expected inverse temperature of the system. 

Finally, we stress that the possible coupling with thermostats set at a temperature 
T can be realized straightforwardly by selecting a subset of links r' C T (or, 
analogously a subset of nodes) and by extracting the pertaining kinetic energies 
according to the Boltzmann distribution, being (3 = 1 /T [3T] . 

4. Metrization 

As foresaid, an interesting characterization of the Ising model on the SG may be 
accomplished by a configurational analysis defined in the wider context of Cellular 
Automata on graphs. Such an analysis can be realized by a particular metrization 
referring non directly to the configuration space, whose points are the states a, b, c, .. 
of the system, but to a peculiar partition space containing - among other elements 
- the cluster distributions of the system. The mathematical framework, which is 
summarized in the Appendix A.l, requires that the graph is endowed with the 
structure of a probability space. Precisely, we consider the triple (G, Q, [i), where the 
measure /i on the subset algebra Q of the gasket G is simply given by the normalized 
number of nodes in every subset. A configuration (or state) on G is a function assigning 
to each node a value in an alphabet K. The set of all possible configurations will be 
denoted as C(G), the configuration space. Since we consider the Ising model on G, the 
alphabet is binary, but all we are going to say is independent of the number | K | . The 
adjacency matrix, combined with the list of values on nodes, allows an easy definition 
of clusters on an arbitrary structure (the procedure is identical to the recognization 
of connected subset in graph colouring): two homogeneous nodes belong to the same 
cluster if they are connected through a path of nodes sharing the same value. Thus, 
every state a determines in a natural way a partition of G, i.e. an exhaustive collection 
a = $(a) = {Ai, A 2 , A n } of disjoint subsets A i} each connected and homogeneous, 
commonly called atoms of the partition. The set of all partitions of G constitutes 
the partition space Z = Z(G). The application $ : C(G) — > Z(G) is many-to- 
one, since, for instance, permutations in K produce different states but the same 
partition. Obviously, Z(G) contains much more partitions than those derived from 
clusterization, e.g. since, in general, atoms do not require to be connected sets, as 
clusters are. In the present case, Z(G) is discrete and finite because G is such, but 
the formalism applies in abstract probability spaces (see [25l l26l l27l l28l [29] ). 

Clearly, when a dynamics is defined on the graph, this determines configuration 
orbits {a(t n )} starting from any initial state a(io), and the corresponding partition 
orbits {a(t n )}, where a(t n ) = $(a(i n )) 

Basic operations between two arbitrary partitions a and (3 are the minimal 
common multiple 7 = a V /3, and the maximal common factor (m.c.f.) or intersection 
a = a A (3 (see Fig. IA1I and Appendix AA for details). The entropy H(a) of the 
intersection a, when calculated between partitions at next steps along an orbit, is an 
index of the relevance of the (instantaneous) non evolving part. 

The metrization of the partition spaces is based on the Rohlin distance, which 
describes the non-similarity of two arbitrary partitions a and f3. This distance, 
requiring the Shannon conditional entropy H(a\j3) of measurable partitions (see 
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Eq. IA.2I) , is given by the functional of Eq. IA.31 we anticipate here: 

d R (a,P) = H(a\P)+H(p\a). (3) 

In order to amplify non similarity, in the Appendix A.2 we present also a method, 
referred to as "reduction process" and denoted n, which acts on couples of partitions 
and uses both operations V and A. More precisely, given two partitions, say a and 
/3, their reduction is obtained by first defining their intersection a = a A (3 and by 
keeping from both partitions only those subfactors a/~ and prime with a, namely 
such that Sfc A a = /3k A a = v. Then, the reduced partitions are given by, a — VjCtj 
and (3 = Vjfij respectively (see Fig. IA21 lower panel). The process ir gives evidence of 
the essentially different sub-partitions of any couple in Z(G), and therefore amplifies 
their distance: du(a, f$) < dn{ir{a, (3)). Hence, by comparing the distance between 
reduced and non-reduced couples, it is possible to introduce an amplification ratio, 
that is dn(Tr[a, f3])/dn(a, /?), which provides further information about the cluster 
distribution and mobility. However, all this analysis has to be performed in correlation 
with other observables. As explained in Appendix A3, the reduction process is 
effective, namely gives rise to a large amplification ratio, whenever one of the two 
partitions, say a, displays at least one cluster which in j3 is exactly decomposed 
into smaller ones. Now, in the case under study a and (3 are partitions defined by 
cluster configurations at two successive steps, and the existence of a common cluster 
is furthered by the special topology of the SG: clusters corresponding to subgraphs 
which are (combinations of) gaskets of generation < g are rather stable (the border 
is made by two vertices only), nonetheless internal fluctuations may occur and hence 
decompose the cluster itself. 

It is worth underlining that while the 7 = a V /3 operation is rather trivial, the 
a = a A f3 operation and the reduction process 7r are quite tricky. Reduction in 
particular is the main algorithmical obstacle in handling large graphs (g > 6). 

Of course, the Rohlin distance is deeply different from the well known Hamming 
distance du in the configuration space, i.e. 

1 N 

dff (a, b) = — p((H,h), (4) 
1 

where a, and bi are the values of the i— th node, p is a distance functional in the 
alphabet. The simplest distance in K is p(x,y) = 1 — 8 xy , leading to Hamming 
distance 0.5 for purely random configurations. Looking for instance to Fig. IA2I (upper 
panel), we would get the maximal Hamming distance for the configurations a and 
b, and a minor distance for a and c, while the Rohlin metrics on the corresponding 
partitions gives null distance in the former case and a high distance in the latter. In 
other terms, d R and dn are deeply different in principle not only because they refer 
to different objects (defined in Z(G) and C(G), respectively), but also because the 
former has do to with mutual distribution of clusters, which could involve geometrical 
features and long range correlations, while the latter is the sum of strictly local 
differences. Therefore, even if the correspondence $ : C — > Z is many-to-one, the 
loss of information should be compensated by the fact that by dn we get a global 
estimate on cluster distributions, instead of a bare sum of uncorrelated differences. 
Finally, more details concerning the implementation of the algorithm for measure of 
the above mentioned metric observables can be found in Appendix AA. 



Metric characterization of cluster dynamics on the Sierpinski gasket 



8 



5. Numerical Experiments and Results 

We now focus on the Ising model on the SG, with the microcanonical dynamics 
described in Sec. [3] The equilibrium regime is ensured by coupling the system to 
thermostats set at the same temperature; in this particular case the simplest way 
is to consider, for a given generation g, the external triangle, i.e the 3 x 2 9_1 links 
defining the perimeter of the gasket; if then we exclude the six angular links, we get 
3 separated thermostats. This could be useful in the future, as it consistutes a quite 
simple way of imposing temperature gradients to the system. However, it should 
be noticed that, due to the fact that the contact between the thermostats and the 
system gets vanishingly small with respect to the bulk as g is increased, the time for 
thermalization is expected to grow with the size of the system. Anyway, as mentioned, 
we are now interested only in the equilibrium behavior. 

Before proceeding, we underline that our checks strongly confirm that the 
thermalized states reached by the microcanonical dynamics are consistent with those 
expected from a canonical dynamics, e.g. based on the Metropolis algorithm. In 
particular, we verified that, for a given temperature, macroscopic observables like 
the magnetization and the energy measured with the two kinds of dynamics are 
indistinguishable . 

For any given temperature T, the geometric observables have been calculated 
as time series starting after a thermalization time, and lasting an observation time 
£max, where time is measured in units of Monte Carlo (MC) steps; i.e. N elementary 
moves. In particular, we consider finite segments of trajectories {a(ifc)} in C or 
{a(tk) = $(a(ifc))} in Z, for tk — 0, 1, 2, ...,t max , as well as the related intersections 
cr(ifc) = a(ifc) A a(tk+i) and the couples of reduced partitions at successive time steps 
(a, 0) — ir(a, /?). From such equilibrium trajectories we obtain segments of time series 
for the following quantities: 

(i) the entropy H(-), applied either to the orbits or {c(ife)}; 

(ii) the Rohlin distance duitk) = dii(a(tk-i),ct{tk)), (see Eq. IA.3|k 

(iii) the amplified Rohlin distance dn(tk) = dn(ir[a(tk-i), a(tk)]); 

(iv) the amplification ratio dn(tk)/dii(tk)', 

(v) the Hamming distance dn(tk) — df/(a(£fc_i), a^,)), (see Eqs. Bland IA.4() . 

Here we focus to the case of unitary time steps (tk+i — tk equals one MC step), 
leaving the study of the role of the time gap to future investigations. 

After setting the experimental parameters (size, thermalization time, confidence 
length of trajectories i ma x), we calculate time averages, variances etc. for each of 
the series above. Results found for different choices of size (we especially focused 
on gaskets of generation g — 4, g = 5 and g = 6, corresponding to N = 120, 
N = 366 and N = 1095 sites respectively) and of t max are qualitatively in very good 
agreement. Moreover, to approach the thermalized state, the initial configuration 
is taken ferromagnetic; this minimizes the likelihood of pinning effects during the 
evolution 30J. In the following we report only the essential information, dropping 
redundant numerical outputs. 

In general, the observables analyzed highlight the existence of different regimes, 
demarcated by remarkable temperatures Xj sa 1.3 and Tjj « 2.0. More precisely, Tj 
corresponds to the flex in the Rohlin distance and to the peak in the amplification 
ratio, while Tu corresponds to the peak in the variance of Rohlin distance and to the 



Metric characterization of cluster dynamics on the Sierpinski gasket 9 




Figure 2. (Color on line) Mean (left panel) and variance (right panel) for Rohlin's 
distance cIr as a function of the temperature; three different sizes are compared. 




Figure 3. (Color on line) Mean (left panel) and variance (right panel) for 
Hamming's distance du as a function of the temperature; three different sizes 
are compared. 

crossover in the intersection entropy (see Figs. 2 — 4). Interestingly, Tj and Tjj recover 
the "critical" temperatures T± and Ti evidenced by thermodynamic analysis (see 
Sec. 2 and [TJH1H]). Indeed, consistently with thermodynamic results, the highlighted 
regimes correspond to a long-range order region and to a disordered region with a 
critical-like transition region in between. More precisely: 

• for T < Tj, (to 2 ) is close to 1 (see [8]), while the distances du and da are 
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approximately (see Figs. 2 and 3); in fact, at such small temperatures a 
ferromagnetic order is established over large length-scales: clusters are constituted 
by few single spots in the large "sea" of equally oriented spins and spin-flips are 
rather unlikely to happen. So, configurations - and partitions - at consecutive 
steps differ for such small spots that distances are extremely small. Moreover, 
during a MC step, the rare spin flips occurring are yet able to change the 
atoms in such a way that overlaps between consecutive partitions are quite 
improbable. This inhibits the reduction, and the amplification is close to 1. For 
T approaching Tj, such spots get larger, but, due to the slowness of the evolution 
at low temperature, their borders can remain sufficiently unchanged for several 
steps. Hence, the emergence of spots within such clusters allows the reduction, 
processing couples of consecutive partitions, to get more effective, and the time 
average of the amplification is manifestly increased. The growth of the spots can 
be retrieved by the intersection entropy increase, as shown in the inset of Fig. 4 
(left panel). 

• for T m Tj, (m 2 ) is close to 0.5 and both djj and d^ start to be significantly 
larger than 0. This is a consequence of the fact that spin flips are getting more 
frequent. More interestingly, the amplification ratio reaches a maximum. We 
have seen that, at the middle of the previous regime, two conditions cooperated 
to start the growth of the amplification ratio: non empty intersection er, and 
the fact that at consecutive steps there exist uniform large clusters which are 
decomposed internally, yielding an effective reduction (as explained in Appendix 
A3). Here, the peak in the amplification ratio means that cluster sizes and the 
speed of the dynamics optimally fit such conditions; moreover, fixed clusters are 
more likely. The exact determination of the peak temperature is tricky, due to the 
growing complexity of the configuration dynamics, but we argue that Tj (or Tj.) 
is a good approximation. This behavior is consistent with the apparent transition 
occurring at T\ as a long-range order breakdown. 

• for Tj < T < Tjj, (m 2 ) is approaching zero. Clusters get more intricate due to 
the growing temperature, and H{o~), which measures the relevance of overlapping 
between successive configurations, exhibits a rapid growth (see Fig. 4, left panel). 
The coexistence of fragmentation of large clusters, which is a signature of decay for 
long range correlations, and overlapping may be related to "critical slowing down" 
effects [31]. Clearly, fragmented clusters have a higher probability to overlap, but 
are unlikely to include small spots; As a consequence the amplification ratio is 
still larger than 1, though rapidly decreasing. Note that H{a) is rather far from 
saturation (namely, uncorrelated chaos) and this suggests the persistence of a 
local short-length order. 

• for T ~ Tji, the amplification ratio is practically 1: such a breakdown of the 
reduction, as forementioned, has a completely different meaning than in the case 
of very low temperatures, where overlapping was insignificant for the smallness 
of sparse clusters in the "big sea" of dominant magnetic orientation. Here, due 
to the smallness of the a atoms, it is unlikely that a sufficiently large cluster is 
internally decomposed at the next step (see also the Appendix A. 3). 

At this temperature the variance of d^ has a maximum (whose value scales like 
l/\/~N), and this constitutes a signature that, along the trajectory, the fluctuation 
of the distance is particularly important. Indeed, the complexity of clusters shape 
and their mobility can give rise to wide fluctuations of distances in time. Also, 
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Figure 4. (Color on line) Intersection entropy H(cr) (left panel) and amplification 
ratio dji/dn (right panel) as a function of the temperature. For the latter, 
accurate numerical data are presented only for generation g = 5, due to time 
consuming calculations; nonetheless we have checked that for other values of g 
the qualitative behavior is robust. The inset of the left panel shows a magnification 
of the departure from zero of the entropy; notice the scmilogarithmic scale. 

H(a) exhibits a crossover: from Tu, the growth is due only to fragmentation, 
while the "critical slowing down" is over. Both the persistent growth of H(a) 
and the value of dn (neatly below 0.5) indicate that the complete chaos is far 
from being established at this temperature and for these sizes. Moreover, the large 
fluctuations make the system more susceptible to spin-flips, i.e. energy changes, 
and this is consistent with the peak in the specific heat. 

• for T > Tu , the Rohlin variance decreases roughly as ~ 1/N toward an 
asymptotic value corresponding to the uncorrelated chaos, whose onset may be 
recognized by dn — 0.5. We argue that the fragmentation and the mobility of 
clusters stabilize the behavior of the time series for global quantities like distances. 
We note however that the maximum for the <i#-variance occurs at T ~ 3, well 
above the corresponding maximum of the ci^-variance at T ~ 2. We deduce that 
in this temperature interval big clusters of all sizes, typical of the previous fractal 
configurations, have started their fragmentation, so that the global shape diversity 
between subsequent configurations is stabilizing and this explains why the da~ 
variance decreases. However, the fragmentation process is still slow enough that 
the local matching between different or equal spins at next steps is highly unstable 
in time, allowing for the growth of the <i#-variance. Indeed, spins belonging to the 
inner part of a cluster (even if of small size) result from the dynamics to be more 
stable than those belonging to its periphery. For T > 3 the fragmentation is such 
that both local overlaps and medium length correlations tend to stabilize the time 
behavior. There is therefore a temperature interval exhibiting a subtle interplay 
between correlation length and time stability. In other terms, in the way to the 
chaos we recognize two time scales: one in terms of global similitude (shape) of 
clusters along the orbits, which starts to stabilize just at the end of the "critical 
slowing down" (T sa 2); the other in terms of local overlaps, whose stabilization 
requires a higher temperature (T ~ 3). The direct visual inspection could be 
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awkward for such properties, undetectable also from the mean magnetization. 
6. Conclusions and perspectives 

In this work we have performed a quantitative characterization of cluster dynamics for 
the Ising model defined on a finite Sierspinski Gasket, whose thermodynamic behavior 
is known to be non-trivial. The analysis is based on a set of geometric observables, 
such as Hamming and Rohlin distances, and on reduction operations among partitions 
which allow to detect the evolving and the stable components of clusters. 

The phenomenology evidenced by previous thermodynamic analysis is qualita- 
tively and quantitatively confirmed by the present metric approach, which provides in 
addition a geometric characterization of the anomalous behavior of the system. Indeed, 
we highlight first the existence of two "critical" temperatures: Tj corresponding to a 
peak in the amplification ratio, meaning that there an abrupt change in the clusters 
behavior occurs, and T/j, corresponding to a peak in the Rohlin variance and other 
crossovers, evidencing the loss of short range order. More precisely, the amplification 
ratio gives detailed information about the development of clusters, starting from lit- 
tle independent spots (amplification ratio equal to 1) to first overlaps (amplification 
ratio larger than 1 and growing to the maximum) while the next phase of decreasing 
ratio indicates the complex effect of decreasing efficiency of the reduction due to the 
fragmentation, up to Tu when the fragmentation is such that, notewithstanding the 
large intersection, the reduction process is inhibited by the extreme improbability of 
clusters with a decomposable inner part (see |Appendix A.3[ ) . 

Another new information we obtain in the high temperature regime is the 
existence of different scales in the destruction of local and correlated order, as 
evidenced by the Hamming and Rohlin variances, with an anticipated peak for the 
latter quantity. 

Thus, the analysis above enabled us to distinguish different behaviors and phases 
demarcated by Ti and Tu as main milestones; we cannot speak of course of "phase 
transitions" , since the onset of these distinct behaviors seems to be a rather smooth 
process. However, from such a detailed description of the collective motion we get an 
insight on the interplay between dynamic and statistical features, especially the decay 
of spatial and temporal correlations. 

This work opens a lot of possible extensions and further insights. The next step 
is the metric characterization of clusters dynamics in condition of non-equilibrium; a 
kind of analysis allowed by the chosen dynamics even in the presence of a disordered 
coupling pattern [21]. 

In particular, it could be interesting to deepen the effect of the double scale 
in the approach to uncorrelated chaos on the conductivity. In the same context, 
another interesting item to explore is the relevance of topology for the cluster diffusion. 
Indeed, the whole set of operations performed on the SG, from cluster identification 
to reduction, immediately applies to any automaton on arbitrary, connected graphs: 
the process depends only on the adjacency matrix, on a generic alphabet and 
on a dynamics, working as a proper external engine generating a succession of 
configurations. Therefore, by the bare substitution of the adjacency matrix, the 
algorithm is ready to fit a great variety of statistical models (e.g. Pott's model 
on lattices of arbitrary dimension or graphs), or other network problems where a 
dynamics is defined. If necessary, nodes could also be weighted, defining alternative 
probability measures. True algorithmic problems and non trivial extensions could only 
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arise from alternative definitions of partitions (assuming e.g. atoms of a different kind 
with respect to the clusters) or from a different factorization, modifying the reduction 
process. 

Appendix A. Technical remarks on metrization 

Appendix A.l. Generalities 

The formalism and general results for partition spaces and Rohlin metrics may be 
recovered e.g. in [25l [26j [27l [28j [29] . Let (M, M,n) be a probability space, that is an 
arbitrary set M, a cr-algebra A4 of subsets of M, and a normalized measure fi on M.. 
In our case the set M is just given by G. 

A partition of M is a finite collection a = (A\, A2, A n ) of measurable disjoint 
subsets covering M, i.e. A t fl Ak = and UkAk = M. The {Afc}'s are called the 
"atoms" of a. The set of all finite measurable partitions is denoted Z = Z(M). The 
unit partition v consists of the single atom M. A partial order in iJ, i.e. a relation 
a < ft, means that ft is a refinement of a; equivalently, every Ak is exactly composed 
with some Bj included in ft. In such case, a is said to be a "factor" of ft. Clearly, 
v < a for every a. 

Such terms as "unit" and "factor" depend on a commutative and associative 
pseudo-product, or composition, 7 = a V ft, denoting the less refined of all partitions 
greater or equal to both a and ft, whose atoms are the non empty intersections of 
the a and ft atoms. If not ambiguous, we can also write 7 = aft. Obviously, ar/ = a 
whenever 77 < a, and in particular av = a for every a. Such properties make the 
result of this operation a kind of "minimal common multiple" . 

Conversely, a — a A ft is the greatest partition such that a < a and a < ft. In this 
case, al\v = v for every a, and ahft — v implies that a and ft are "relatively prime" 
(i.e. they have no common factor) . Therefore, the result is a sort of "maximal common 
factor". See Fig. I All for an example of product and intersection among partitions. 

A partition may represent a probabilistic experiment with non overlapping 
outcomes A%,..,An, where the "atomic" event Ak has probability fj,(Ak). A factor 
is therefore a sub-experiment of the finer experiment, grouping several outcomes 
as equivalent: for instance, "odd or even" is a two-atoms sub-experiment of the 
{1, 2, 3, 4, 5, 6} dice experiment. 

The Shannon's Entropy H{a) defined on every partition as 

n 

H{a) =-J2»( A i) ln K^i) , (A.l) 

i=l 

is a measure of the mean information obtained from the experiment. If j3 = 
(£?i, B m ) is another partition, the conditional entropy of a with respect to ft is 

H(a\ft) = - E E ^ n Bk) ^ M 7^7 j , (A.2) 

i=l k=l k > 

where, as usual, one takes xh\x = for £ = 0. This conditional entropy is the 
mean residual information obtained from a when the result of ft is known. Note that 
the Shannon entropy depends only on the distribution of the atom measures, not on 
their nature or "shape" (this term coud have no meaning in abstract spaces). On 
the contrary, the mutual relations among atoms (and possibly their shapes) directly 
influence the conditional entropy (see Figure A.2, upper row). 
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Figure Al. Product (upper panel) and intersection (lower panel) between two 
examples of partition realized on a square lattice. 



Now, the Rohlin distance dji is defined in Z(M) by 

d R (a,P)=H(a\l3)+H(P\a), (A.3) 

and it may be considered as a measure of the overall non-similarity between a and j3. 

If M is finite, a configuration or state a on M is a function assigning to each 
point ijEMa value = f(xi) in an alphabet IK. All possible configurations form a 
space C = C(M). In C the Hamming distance du is defined by 

d H (a,b)=AfJ2p(^h), (A.4) 

i 

where p{ai, hi) is a distance in K and Af a possible normalization coefficient. 

To each configuration corresponds an exhaustive collection of clusters, i.e. 
connected subsets of M with homogeneous value in IK, defining a particular partition in 
Z(M). This establishes a many-to-one correspondence $ : C — >• Z, making possible the 
comparison between ci//(a, b) in C and d R (a, f3) in Z, where a = 3>(a) and (3 = <I > (b). 



Appendix A. 2. Reduction 

The essential non-similarity between two partitions could be confused and weakened by 
the presence of a tight common factor, i.e. a common sub-partition (see Fig. lA2[ lower 
panel). Therefore, we would eliminate common factors as far as possible, a "reduction" 
which is expected to amplify the distance. However, this operation (analogous to 
the reduction to minimal terms for fractions [17]) is not uniquely definite because 
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Figure A2. Upper row: The three configurations depicted evidence the difference 
between Hamming and Rohlin distances: du is maximal for configurations a and 
b, and minor for a and c; dft is null in the former case and large in the latter. 
Lower row: For the partitions a and ft the common factor is given by <r, while 3 
and P are the reduced partitions respectively 



partitions, differently from integers, do not admit a unique factorization into primes. 
The role of primes (i.e. indecomposable) factors can be played by dichotomic sub- 
partitions, which are still extremely redundant (2 n_1 — 1 indeed for a partition with 
n atoms). 

For a = (Ai, A2, A n ) we shall define therefore a restricted family E(a) of 
"elementary" dichotomic factors a\, a.2, S n such that 

(i) E(a) must be well defined for every ae2; 

(ii) E(a) does not contain more than n (the number of atoms in a) elementary factors; 

(iii) v fc=i"fe = a - 

A universal choice consists in taking as dichotomic factors 5^ = (AktA^.), the 
partitions formed by single atoms and their complements to M. Elementary factors 
of this form, used throughout in the present paper, will be called "simple" . 

For a couple a and (3, once their elementary factors E(a) and E(/3) have been 
defined, the reduction process consists in the following steps: 

(i) define the maximal common factor a = a A ft; 

(ii) drop from E(a) and E(/3) those factors which are not relatively prime with tr, and 
note the surviving factors ctk and j3j respectively (this means SfcAcr = j3j Act = ii); 

(iii) define a = VfcSfc and (3 = Vjjlj. 

The reason of step ii, which seems to be cumbersome with respect to the simple 
dropping of common factors in E(a) and E(/3), is that in general two families could 
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Figure A3. Upper row: a and /3 are two partitions of the square, each of four 
atoms; a is their intersection or m.c.f. of their atoms. Lower row: List of relevant 
elementary dichotomic factors, where the black specifies the atom and the white 
the complementary set. 



a 



a 



Figure A4. The partitions & and are the reduced partitions of those appearing 
in figure [A3] atoms are now three, individuated by black, white and grey. The 
partition a is the reduced intersection. Notice that the grey atom is non connected. 



have no common factors and, nevertheless, a ^ v. This happens, for instance, when 
a < j3 with no common elementary factors. Then a — a and a = v with the reduction 
above, while a = a with the dropping of common factors. 

It results that dn(a, (3) > dn(a,/3), as requested [T71 [TBI [25] , 

Note however that while two relatively prime partitions are already reduced, not 
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necessarily two reduced partitions are relatively prime: see for instance the case of 
Figure A. 3, which could represent small portions of wider partitions deduced from 
almost chaotic configurations. In the next subsection we give more precise details on 
this. 

The correspondence it : (a, (3) — > (2,/?) is many-to-one and idempotent, i.e. 
7r o 7r = 7T. It is a sort of projection from Z x Z on the subset of irreducible pairs. 

The reduction process, therefore, essentially depends on the choice of the family 
E(a) of elementary factors. Besides simple factors, other families exist, which 
in particular cases could prove more convenient for algorithmic reasons or for the 
observer's attitude in the probabilistic experiment. While simple dichotomic factors 
correspond to looking at the "occurrence or not" of single atomic outcomes, other 
attitudes could isolate dichotomic outcomes enjoying supplementary properties which 
are not universal but depend, typically, on some additional geometrical structure of M 
(order, connection, orientation etc.). For instance, in previous works on rectangular 
lattices jTHl [20] elementary factors were identified by external contours of clusters, 
a choice intended to optimize the simple connection of the factor atoms. This is 
not convenient on general graphs, where M = G, because the determination of 
external surfaces could be cumbersome or impossible. Therefore, in the present work, 
elementary factors will always be the simple factors {Ai, Af\. 

Appendix A. 3. Reduction vs. Amplification 

Let a and (3 such that a — a A (3 ^ z/, and a — (Si, S2, S q ). Every Sk € cr is the 
union of some subsets {Aki} and {Bkj} of atoms in a and j3. Possibly, such subsets 
may be of one single atom or more: we indicate s£ and m£ the single or multiple 
atoms cases for a, and analogously s^ and m? k for (3. Clearly, at fixed k, s% = Ak 
and = Bk, but it is useful to keep a distinct notation in order to remember that 
such atoms are not only in the partitions but also in their m.c.f. Therefore, Sk may 
be composed in four forms: (s^,s^), (s^,m^), (m^,s k ) and (m^,m^). For instance, 
in Figure IA"2l the m.c.f. a has one (m°,s^), one (s a 1 m^) and three (s Q ,s' 3 ) atoms, 
while, in Figure IA"3"1 a has one (m a ,m 13 ) and two (s Q ,s^) atoms. 

Proposition 1: the atoms of a are all those of a contained in the {m^} groups, plus 
one atom constituted by the intersection of their complementary sets, or equivalently 
by Ufcs£ (and similarly for (3). 

The proof immediately follows from the fact that the (Ak, A^) elementary factors 
are dropped in the reduction process if and only if Ak = sj?. This is true in any abstract 
partition space, whenever simple factors are used; notice however that, speaking of 
partitions generated by connected configurations, as those considered in the present 
paper, this last atom could be non-connected (an example in the figure |A"4| . 

Proposition 2: the atoms of a = a A j3 are: i) all the Sk G cr of the (mj?,mf) form; 
ii) one more atom (the complementary part to the union of the previous ones) if at 
least one term (s£,s^) exists; or in) two more atoms if, besides the (m^,m^), only 
mixed terms (s£,mf) and (?n£,sf) exist. 

The point i) follows from Proposition 1, since the {m^} and {m^,} groups 
reconstitute common subpartitions in a and /3; point in) depends on the fact that, in 
absence of (s£,sf) terms, there is an exact correspondence between U^s^ and U^m^! 
(and equivalently UfeS^ and U^m^ ), where the k index runs over the mixed-type 
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Sfc's only; the two supplemetary atoms are therefore U^s^ and U^s^; as to point ii), 
we observe that the complementary set of the {m%}, i.e. U^s^, cannot be exactly 
decomposed by the atoms of {m^,}, and vice versa, just because the presence of the 
supplementary simple terms. 

Several easy corollaries follow. For instance, when simple factors are used, 
Propositions 1 and 2 constitute a direct constructive proof that the reduction tt is 
a projection. Or that there is no reduction at all when in a there are only (m^,m^) 
terms, or when a has only two atoms. Or else that in order to have a = v 1 no terms 
(m£, mt ) are allowed and at least one term (s£, ) is necessary. Moreover, if for one 
partition, say a, no {m^} group exists, then a — a and a = v. Etc. 

An important item is the relation between reduction (meant as pattern simpli- 
fication, as stated in Proposition 1) and the amplification ratio dn{a, /3)/dji(a, /?), 
measuring the metric effectiveness of the reduction. 

Proposition 3: The necessary and sufficient conditions in order to have dn(a, j3) > 
dn(a,f3), and therefore amplification ratio > 1, are the existence in a of at least one 
atom of mixed form, (sj*,m^) or (m^,s^), and the existence of at least two distinct 
group s% from the same partition (x = a or /3). 

This may be proved using the important equality 

d R (a,/3) = 2H(aV/3)-H(a)-H(/3), (A.5) 

(see [17l [25] ) . In the computation of all entropies appearing in Eq. IA.51 one can split 
the Shannon's sums (See Eq. IA.1|) along the atoms Si, S2, S q of a, because the 
intersections AiC\Bj appearing in aV/3 are surely empty for different Sk's. Therefore, 
using Propositions 1, for atoms of the (mj?,m^) form, the partial contributes to 
entropies in Eq. IA.5I are the same in dn(a,f3) and dn(tt,p). In absence of mixed 
terms, the remaining part of both a and ft is constituted by the same atoms = st 
, and the contribute to the distance is 0, exactly as for the supplementary single atom 
UfcS^ = UfcS^ common to a and (3. 

Assume now that also mixed terms (s^,m?) and (tj^,s^) are present (and 
possibly also terms (s%, sf )) in such a way that at least two s% or two s^ appear. Using 
again Eq. IA.5I for these components, H(a V (3) and H(a V j3) are equal, because the 
intersection of the multiple atoms in one partition give the same result with the single 
atoms in the other partition (before reduction) or with their union (after reduction). 

As to the subtracted quantities in Eq. IA.51 the contributes from H(a) and H((3) 
are strictly greater than H(a) and H(f3) because there are contributes from separate 
atoms in the former case, and from their union in the latter (this is why at least two 
s£ or two sj? are requested). Clearly, elementary convexity properties of the — x\nx 
function used here, see e.g. [2"5] , Therefore dn(a,ft) is strictly greater than d fl (a,/3), 
and after the subtraction in Eq. IA.51 the distance is increased, i.e. the amplification 
ratio is greater than one. 

Proposition 3 clarifies that the metric effectiveness of the reduction, besides 
the pattern simplification, is due to the difference between — J2k A*( S I) m M s fe) an d 
— /i(UfcS^) ln/i(UfeS^), where 1 is a or /3. In particular, if a and /3 are partitions 
defined by cluster configurations at two successive steps as those considered in the 
present paper, then the amplification requires that there exist big clusters at one time 
exactly decomposed at the next (or previous) time into smaller ones. 
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Appendix B. The algorithm 

The N x N adjacency matrix A of the graph G and its state a = (01,02, ...,ajv), 
i.e. the K- valued list of sites labelled (1,2,..., iV), constitute the essential information 
necessary to work out the partition algebra and the metrization algorithm, which is 
independent of the way the states are generated (dynamics). 

By standard "colouring" techniques, connected sites are iteratively recognized 
testing their value in K, and the partition is easily produced as a list of grouped 
labels, one group for every atom. For a ten sites graph, for instance, a particular 
partition could be a = {Ai, A 2 , A 3 } = {(1, 3, 4), (2, 10), (5, 6, 7, 8, 9)}, with measures 
3/10, 2/10, 5/10 for ^1,^2,^3 respectively. A useful representation for the atoms 
is a binary string of length N: e.g., in the example above, A\ = (1011000000), etc. 
Note that from now on the topological nature of the graph does not influence the 
operations, which regards only the label lists. The "simple factors" of the family E(a) 
are immediately defined. In the example, a\ = {(1, 3, 4), (2, 5, 6, 7, 8, 9, 10)}, etc. (the 
corresponding binary strings are obviously complementary). Clearly, the extreme 
simplicity of this procedure could be replicated with other choices of dichotomic 
factors, e.g. by taking the "internal-close" and the "external-open" parts of Jordan 
contours (provided that such contours are well defined). Only in the very special case 
of one-dimensional chains the latter choice (internal-external, i.e. left-right) proves not 
only more effective but easier than the former, because contours are left extremities 
of semi-open intervals. 

The V operation may be easily implemented by boolean intersections on the atom 
lists. This is enough for the Rohlin distance computation by formula lA.ll thanks to 
the equality of Eq. IA. 51 above. 

More attention is required for the A operation, which is the key step to define the 
maximal common factor a = a A/3, but the task may equally accomplished by boolean 
operators: every atom of a is built site by site, testing the simultaneous belonging 
to some atoms of a and /3, up to exhaustion. Once the maximal common factor is 
defined, most of the computation time is spent in the storage and management of the 
surviving simple factors a,k and j3j, along the criteria of the subsection | Appendix A.3 
However, all this does not imply conceptually new operations. Reduced partitions, 
and their atoms, which in principle are defined by intersections on the atoms of such 
surviving factors, can take advantage of Proposition 1 in subsection |Appendix A.3[ 
completing the ir process. For the amplified distance da, the functional in Eq. IA.5I is 
applied to the couple (otk,Pj), and the observations of the previous subsection reveal 
useful again. 
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